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The cumulative probability distribution of sparseness time interval in the Internet is 
studied by the method of data analysis. Round-trip time between a local host and a 
destination host through ten odd routers is measured using the Ping Command, i.e., 
doing echo experiment. The data are found to be well described by the ^-exponential 
distributions, which maximize the Tsallis entropy indexed by q less or larger than unity, 
showing a scale-invariant feature of the system. The network is observed to itinerate 
over a series of the nonequilibrium stationary states characterized by Tsallis statistics. 
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The Internet is a complex system, which has highly intricate tangle, cluster and 
hierarchical structures, strong spatio-temporal correlation with feedback, self- 
organization, and connection diversity. The structure emerging from the actions of a 
large number of the users may efficiently be understood within the statistical 
mechanical framework and its suitable generalizations. For example, a work in [1] 
reports emergence of scaling behavior and the associated power-law distribution of 
connectivity of nodes. These concepts are known to be essential for the network to be 
resilient and robust to random errors, breakdown, and attack [2-4]. 

From the statistical and dynamical viewpoints of the network, of particular interest 
are the stationary states under nonequilibrium conditions. Tsallis statistics [5] based on 
a nonextensive entropy [6] aims to offer a theoretical basis for analyzing complex 
systems at such states. It has successfully been applied to a variety of problems 
including anomalous diffusion [7,8], Levy flight [9-11], fractal random walk [12], 
complex high-energy processes [13-17], cosmic rays [18], turbulence [19], earthquakes 
[20], stock markets and incomes [21,22], nonlinear maps at the edge of chaos [23-29], 
stochastic resonance [30], protein folding and biomolecules [31,32], citation networks 
of scientific papers [33], urban agglomeration [34], and linguistics [35]. 

In this Letter, we present an experimental evidence that Tsallis statistics in fact 
describes the scale-invariant stationary states of the Internet. 

"Echo experiment" we have performed uses the Ping Command [36,37]. A Ping 
signal is emitted from a local host computer, takes a round trip to a destination host (i.e., 
a site accessed), and returns to the local host through ten odd routers. The route of the 
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signal emitted to the destination is fixed and traced. Each router is connected with the 
whole network in a time-dependent manner. The next signal is sent immediately after 
the previous one returns. Such a time interval is typically < 130 ms and is not included 
in our data analysis. Using the whole collected data of the echo experiment, the 
threshold value indicating congestion is appropriately defined (see the later discussion). 
Actually, the result turned out to be not sensitive to the definition of the threshold. We 
have calculated each time interval given by the amount of the round-trip time below the 
threshold value between two successive thresholds. This interval is referred to here as 
"sparseness time interval" denoted by T. Observation shows that the Internet itinerates 
over a series of the stationary states, which are all described by the Tsallis ^-exponential 
cumulative probability distribution of T, indicating a scale-invariant feature of the 
system. In particular, both the < q < 1 and q > 1 cases occur and nothing is special in 
the limit q — > 1. Regarding sparseness time interval rescaled by its average, the more 
congested the network is, the smaller value q takes. In this sense, the entropic index 
characterizes the degree of congestion. 

Before presenting the experimental results, let us briefly summarize the basics of 
Tsallis' scale-invariant statistics. This theory aims to offer a framework for describing 
statistical properties of complex systems at their stationary states based on principle of 
maximum nonextensive entropy. In the present case, the fundamental random variable 
is sparseness time interval, T. p(t) dx is the probability of finding the value of 
sparseness time interval in the range, [t, t + dx\ . Then, p(r) is a stationary probability 
distribution in Tsallis statistics if it optimizes the Tsallis entropy [6,38] 
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(i) 



under the constraints on normalization 



Jjr p(r) = 1 



(2) 



and the normalized ^-expectation value [39] of sparseness time interval 



<?> q =\dTT P q (t). 



(3) 



Here, q and o in Eq. (1) are the positive entropic index and a scale factor of the 
dimension of time, respectively. P q {t) in Eq. (3) is the escort distribution [40] defined 
by 



p"{z) 



\dt' p q (T')' 



(4) 



The optimal distribution is calculated to be 



p(*) = ^e/--(*-<T>J 



z q 'V c 



(5) 
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T max / O 

Z q = J" dtej 9 

o 



V c 



(6) 



c= )^-[ap{T)Y 
J a 



(7) 



f3 in Eqs. (5) and (6) is the Lagrange multiplier associated with the constraint in Eq. (3). 
e q (x) stands for the "^-exponential function" defined by 



[l + (l-q)x\ m ~ q) (l + (l-q)x>0) 
(l + (l-q)x<0) 



(8) 



whose inverse is the "^-logarithmic function" 



x l ~ q -1 

ln 9 W = ^ 

1 -a 



(9) 



Accordingly, T max — > °° if g>l, whereas T max = T /(l- q) if 0<g<l. Here, 
T =[c + (l-g)/3<T> 9 ]//3, which can be shown to be always positive [11]. p(r) is 
recast into the following form: 



P(*) = 



(10) 
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p(r) is seen to be the Zipf-Mandelbrot distribution with a heavy tail if q > 1. Both the 
normalizability condition and finiteness of < T > in Eq. (3) require the entropic index 
to satisfy g<2. 

In the limit g — > 1, the Tsallis entropy converges to the Boltzmann-Shannon entropy, 

S = dr p(r)\n[a p(r)], and correspondingly p(r) becomes the Boltzmann-type 

exponential distribution since, in this limit, e q (x) and \n q (x) approach to the ordinary 
exponential and logarithmic functions, respectively. However, such a limit does not 
play any special roles in the present work. 

An important point in Tsallis statistics is that the quantity to be compared with the 
observed distribution is not p(r) in Eq. (10) itself but its associated escort distribution 
[41]. Therefore, the cumulative probability distribution should be defined by 

P( > T ) = JJ- dx' P q (f). From Eq. (10), it is found to be given by 

P(>t) = «,(-t/t ). (11) 

Below, we discuss how the cumulative probability distributions of this form are realized 
in the Internet. 

In Fig. 1, we present an example of an observed time series of sparseness time 
interval, % . Three distinct stationary regimes, a, b, and c, may be recognized. (Strictly 
speaking, identification of stationary state depends on time scale. Here, we are 
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employing user's typical time scale, i.e., 10 minutes ~ 1 hour.) In Fig. 2 a, b, and c, the 
corresponding cumulative probability distributions of t are plotted on the log-log scale. 
The threshold value indicating congestion is defined here by the mean value plus half of 
the standard deviation. The experimental data are represented by the dots, whereas the 
curves depict the ^-exponential functions. In particular, the lower ones are drawn on the 
"semi-g-log" scale with different values of q. The resultant straight lines imply that the 
observed cumulative probability distributions are in fact the Tsallis g-exponential 
distributions. 

For comparison, we present Fig. 3 to show that there also exist stationary states, at 
which the values of the entropic index are less than unity. 

These results imply that the network undergoes a series of transition from one 
stationary state to another: [q v T J — > [q 2 , T 2 ) — > [q 3 , T 3 ) — > L . Each stationary 
state is scale-invariant and maximizes the Tsallis entropy. The points of transition 
correspond to a catastrophic changes in the time series of round-trip time (not 
sparseness time), e.g., sudden heavy congestion. 

In conclusion, we have found that the Internet itinerates over a series of the scale- 
invariant nonequilibrium stationary states described by Tsallis statistics. We wish to 
emphasize that the time series of sparseness time is highly nonstationary and non- 
Gaussian. This fact makes it difficult to identify stationary regimes by power spectrum 
analysis, in general. The present work indicates usefulness of Tsallis statistics for 
defining stationary states in the time series exhibited by complex systems. 



7 



Acknowledgment 



We would like to thank Dr. V. Latora and Dr. M. Takayasu for useful discussions. 
References 

[1] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 
[2] Y. Tu, Nature 406, 353 (2000). 

[3] R. Albert, H. Jeong, and A.-L. Barabasi, Nature 406, 378 (2000). 
[4] R. Cohen, K. Erez, D. ben-Avraham, and S. Havlin, 

Phys. Rev. Lett. 85, 4626 (2000). 
[5] S. Abe and Y. Okamoto eds., Nonextensive Statistical Mechanics and Its 

Applications (Springer-Verlag, Heidelberg, 2001). 
[6] C. Tsallis, J. Stat. Phys. 52, 479 (1988). 
[7] C. Tsallis and D. J. Bukman, Phys. Rev. E 54, R2197 (1996). 
[8] A. Upadhyaya, J. -P. Rieu, J. A. Glazier, and Y. Sawada, 

PhysicaA293, 549 (2001). 
[9] D. H. Zanette and P. A. Alemany, Phys. Rev. Lett. 75, 366 (1995). 
[10] D. Prato and C. Tsallis, Phys. Rev. E 60, 2398 (1999). 
[11] S. Abe and A. K. Rajagopal, J. Phys. A 33, 8723 (2000). 
[12] A. Robledo, Phys. Rev. Lett. 83, 2289 (1999). 



8 



[13] D. B. Walton and J. Rafelski, Phys. Rev. Lett. 84, 31 (2000). 

[14] W. M. Alberico, A. Lavagno, and P. Quarati, Eur. Phys. J. C 12, 499 (2000). 

[15] O. V. Utyuzh, G. Wilk, and Z. Wlodarczyk, J. Phys. G 26, L39 (2000). 

[16] I. Bediaga, E. M. F. Curado, and J. M. de Miranda, Physica A 286, 156 (2000). 

[17] C. Beck, Physica A 286, 164 (2000). 

[18] C. Tsallis, J. C. Anjos, and E. P. Borges, e-print astro-ph/0203258. 

[19] C. Beck, G. S. Lewis, and H. L. Swinney, Phys. Rev. E 63, 035303 (2001). 

[20] S. Abe and N. Suzuki, e-print cond-mat/0208344. 

[21] L. Borland, Phys. Rev. Lett. 89, 098701 (2002). 

[22] E. P. Borges, e-print cond-mat/0205520. 

[23] U. M. S. Costa, M. L. Lyra, A. R. Plastino, and C. Tsallis, 

Phys. Rev. E 56, 245 (1997). 
[24] M. L. Lyra and C. Tsallis, Phys. Rev. Lett. 80, 53 (1998). 
[25] U. Tirnakli, C. Tsallis, and M. L. Lyra, Eur. Phys. J. B 11, 309 (1999). 
[26] F. A. B. F. de Moura, U. Tirnakli, and M. L. Lyra, 

Phys. Rev. E 62, 6361 (2000). 
[27] V. Latora, M. Baranger, A. Rapisarda, and C. Tsallis, 

Phys. Lett. A 273, 97 (2000). 
[28] R. S. Johal and U. Tirnakli, e-print cond-mat/0201464. 
[29] F. Baldovin and A. Robledo, e-print cond-mat/020537 1 , 

to appear in Phys. Rev. E. 
[30] H. S. Wio and S. Bouzat, Braz. J. Phys. 29, 136 (1999). 



9 



[31] I. Andricioaei and J. E. Straub, Phys. Rev. E 53, R3055 (1996). 

[32] C. Tsallis, G. Bemski, and R. S. Mendes, Phys. Lett. A 257, 93 (1999). 

[33] C. Tsallis and M. P. de Albuquerque, Eur. Phys. J. B 13, 777 (2000). 

[34] L. C. Malacarne, R. S. Mendes, and E. K. Lenzi, 

Phys. Rev. E 65, 017106 (2001). 

[35] M. A. Montemurro, Physica A 300, 567 (2001). 

[36] I. Csabai, J. Phys. A 27, L417 (1994). 

[37] M. Takayasu, H. Takayasu, and T. Sato, Physica A 233, 824 (1996). 

[38] S. Abe, e-print cond-mat/0206078, to appear in Phys. Rev. E. 

[39] C. Tsallis, R. S. Mendes, and A. R. Plastino, Physica A 261, 534 (1998). 

[40] C. Beck and F. Schlogl, Thermodynamics of Chaotic Systems: An Introduction 

(Cambridge University Press, Cambridge, 1993). 

[41] S. Abe and A. K. Rajagopal, J. Phys. A 33, 8733 (2000). 



10 



Figure Captions 

Fig. 1 Time series data of the sparseness time interval taken from 3.50 a.m. (initial 
time) to 6.07 a.m. on 8 February, 2002. From the local host, 
buffalo.matsudo-ap3.dti.ne.jp [203.181.67.200], to the destination host, 
ring.so-net.ne.jp [202.238.95.103], through 11 routers. The curve is drawn 
based on 31675 measured data points. Roughly, three different 
nonequilibrium stationary states, a (3.50 a.m. - 4.15 a.m.), b (4.15 a.m. 
- 5.06 a.m.), and c (5.06 a.m. - 6.07 a.m.), may be recognized. 

Fig. 2 Log-log plots of the cumulative probability distributions associated with the 
the states, a, b, and c. The observed data are represented by the dots, whereas 
the Tsallis distribution by the solid lines. The lower ones are drawn on the 
semi-g-log scale, a: q - 1.07, T = 2.50 x 10 3 ms, and 4373 data points, 
b: g = 1.12, T = 4.35x10 2 ms, and 13587 data points, c: q = 1.16, 
T = 1.00xl0 3 ms, and 13715 data points. 

Fig. 3 An example with the value of q less than unity. The hosts are the same as in 
Fig. 1. Data was taken from 4.42 a.m. to 5.31 a.m. on 13 February, 2002. 
q = 0.73, T = 2.27 x 10 4 ms, and 14897 data points. 
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